SFMAR 


Systems  Center 
PACIFIC 


TECHNICAL  REPORT  2085 
August  2016 


FANS  Simulation  of 
Propeller  Wash  at  Navy  Harbors 
(ESTEP  Project  ER-201031) 


Pei-Fang  Wang 

SSC  Pacific 

Hamn-Ching  Chen 

Texas  A&M  University 


Approved  for  public  release. 


SSC  Pacific 
San  Diego,  CA  92152-5001 


SSC  Pacific 

San  Diego,  California  92152-5001 


K.  J.  Rothenhaus,  CAPT,  USN  C.  A.  Keeney 

Commanding  Officer  Executive  Director 


ADMINISTRATIVE  INFORMATION 

This  report  was  was  developed  in  support  of  Environmental  Security  Technology  Certification 
Program  (ESTEP)  Project  ER-20 1-031  by  the  Environmental  Sciences  Branch  (Code  71750),  of 
the  Advanced  Systems  and  Applied  Sciences  Division  (Code  71700),  Space  and  Naval  Warfare 
Systems  Center  Pacific  (SSC  Pacific),  San  Diego,  CA;  and  the  Ocean  Engineering  Program, 
Zachry  Department  of  Civil  Engineering,  Texas  A&M  University,  College  Station,  TX, 


Released  by 
P.  J.  Earley,  Head 
Environmental  Sciences  Branch 


Under  authority  of 
A.  J.  Ramirez,  Head 
Advanced  Systems  &  Applied 
Sciences  Division 


This  is  a  work  of  the  United  States  Government  and  therefore  is  not  copyrighted.  This  work  may 
be  copied  and  disseminated  without  restriction. 


The  citation  of  trade  names  and  names  of  manufacturers  in  this  publication  is  not  to  construed  as 
official  government  endorsement  or  approval  of  commercial  products  or  services  referenced  herein. 

Dell®  is  a  registered  trademark  of  Dell,  Inc. 

IBM®  is  a  registered  trademark  of  International  Business  Machines  Corporation. 

Intel®  is  a  registered  trademark  of  Intel  Corporation. 

Gridgen®  is  a  registered  trademark  of  Gridgen  Glyph. 

Linux®  is  a  registered  trademark  of  Linus  Tolvalds. 

MATLAB®  is  a  registered  trademark  of  MathWorks. 

TecPlot®  is  a  registered  trademark  of  .Tecplot,  Inc. 


SB 


EXECUTIVE  SUMMARY 


Propeller  wash  induces  high  shear  stresses  on  seafloor  which  may  cause  sediment  resuspension  in 
Department  of  Defense  (DoD)  harbors.  To  improve  understanding  of  the  sediment  erosion,  transport, 
dispersion,  and  re-deposition  processes,  it  is  desirable  that  advanced  computational  fluid  dynamics 
models  should  be  used  to  provide  detailed  resolution  of  the  velocities  and  bottom  shear  stresses 
induced  by  the  propeller  wash  in  confined  shallow  water  basins.  In  this  study,  the  Finite- Analytic 
Navier-Stokes  code  was  employed  to  solve  the  Reynolds-Averaged  Navier-Stokes  equations  in 
conjunction  with  advanced  near-wall  turbulence  model  for  several  propeller-wash  scenarios 
involving  an  Arleigh  Burke-class  destroyer  with  twin-screw  propellers  and  a  tugboat  with  two  ducted 
propellers.  This  model  enables  us  to  evaluate  the  effect  of  water  depth,  ship  speed,  propeller  rotating 
speed,  and  pier  wall  configuration  on  the  propeller-induced  shear  stresses  distributions. 
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1.  INTRODUCTION 


In  this  study,  the  Finite- Analytic  Navier-Stokes  (FANS)  code  was  employed  for  propeller-wash 
study  in  U.S.  Navy  harbors.  The  FANS  code  solves  Reynolds-Averaged  Navier-Stokes  (RANS) 
equations  together  with  advanced  turbulence  models  in  general  curvilinear  coordinate  systems  using 
overset  (chimera)  grids.  The  overset  grid  system  facilitated  the  simulation  of  arbitrary  relative 
motions  among  various  computational  grid  blocks  such  as  those  encountered  wave-current-body 
interactions  and  vortex-induced  vibrations.  The  FANS  code  consists  of  the  following  main 
components:  (1)  finite-analytic  method  for  the  solution  of  compressible  and  incompressible  RANS 
equations  and  energy  equation  in  general  curvilinear  coordinates;  (2)  dynamic  chimera  domain 
decomposition  technique  for  overlapped,  embedded,  or  matched  grids  including  relative  motions; 

(3)  near-wall  Reynolds  stress  (second-moment)  and  two-layer  k- £  turbulence  models  for  turbulent 
boundary  layer  and  wake  flows;  (4)  large  eddy  simulation  for  unsteady  chaotic  eddy  motions, 

(5)  linear  and  nonlinear  wave  effects;  (6)  level-set  method  for  interface-capturing  between  two 
different  fluids;  (7)  detailed  propeller  flow  simulations  or  interactive  coupling  with  propeller 
performance  programs;  (8)  coupling  with  six-degree-of- freedom  motion  program  for  ship,  structure, 
wave,  and  current  interactions;  and  (9)  multi-processor  parallelization  for  large-scale  computational 
fluid  dynamics  (CFD)  applications.  The  combination  of  these  methods  provides  a  unique  capability 
for  modeling  complex  fluid  flow  and  heat  transfer  around  practical  three-dimensional  configurations 
including  viscous  and  violent  free  surface  effects. 

The  FANS  code  has  been  used  extensively  for  the  simulation  of  ship  motions  under  very  shallow 
water  conditions.  Chen  and  Chen  (1998)  and  Chen,  Liu,  Huang,  and  Davis  (2000)  performed 
berthing  simulations  for  DDG-5 1  and  AOE-6  ships  with  underkeel  clearance  ranging  from  5  to  20% 
of  the  ship  draft.  Chen  and  Huang  (2003)  and  Huang  and  Chen  (2003)  also  performed  time-domain 
simulation  of  berthing  operations  involving  a  modular  hybrid  pier  (MHP)  with  2  mooring  dophins, 

2  moored  LHD  ships,  one  berthing  LHD  ship,  a  tug  boat,  7  fenders  and  12  mooring  lines.  In  these 
simulations,  a  very  shallow  water  depth  of  28  ft  was  intentionally  chosen  to  confine  the  underkeel 
clearance  of  the  moored  and  docking  ships  to  1  ft  or  3.7%  of  the  ship  draft  (=  27  ft).  The  FANS  code 
was  also  used  in  Chen,  Lin  and  Huang  (2002a,  2002b)  and  Chen,  Lin,  Liut  and  Huang  (2003)  for  the 
study  of  multiple-ship  interactions  in  dredged  navigational  channels.  More  recently,  Huang  and  Chen 
(2007,  2010)  used  the  FANS  code  for  site-specific  passing  ship  effects  on  a  docked  ship  moored  to  a 
floating  pier  in  Norfolk  harbor.  These  simulation  results  clearly  demonstrated  the  capability  of  the 
FANS  code  to  model  complex  interactions  between  Navy  harbor  facilities  and  its  client  ships  in  a 
real  waterfront  ambience,  including  site-specific  conditions  such  as  sea  bed  bathymetry,  shorelines, 
harbor  geometry,  navigation  channels,  hull  shapes,  and  arbitrary  ship  motions. 

In  addition  to  the  simulation  of  multiple-ship  interactions,  the  FANS  code  was  also  used 
extensively  for  a  wide  range  of  fluid-structure  interaction  problems,  including  propeller-ship 
interactions  (Chen  and  Lee,  2004;  Lee  and  Chen,  2005;  Pontaza,  Chen  and  Lee,  2006),  greenwater 
and  extreme  slamming  of  ship  in  random  waves  (Chen  and  Chen,  2014),  wet-deck  slamming  (Chen 
and  Yu,  2009),  LNG  tank  sloshing  (Chen,  2011),  hurricane  wave  loads  on  offshore  platform  and 
jack-up  structure  (Chen,  2010,  2013),  vortex-induced  vibration  of  deep  water  risers  (Huang,  Chen 
and  Chen,  2010,  2011,  2012)  and  riser  interferences  (Chen,  Chen  and  Huang,  2013),  and  scour 
around  bridges  (Chen,  2002;  Briaud,  Chen,  Li  and  Nurtjahyo,  2004;  and  Briaud  and  Chen,  2006). 
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In  the  FANS  code,  the  unsteady  Reynolds-Averaged  Navier-Stokes  (RANS)  equations  for 
incompressible  flow  in  curvilinear  coordinates  are  formulated  in  general  curvilinear  coordinates 

(£',0  =  (&7,£,0: 


u\  =  0 

+ U %  +  }, + gV,  -  j-egJku:jk  =  o , 


(1) 

(2) 


where  U  and  ii  represent  the  mean  and  fluctuating  velocity  components,  and  g'J  is  conjugate  metric 
tensor,  t  is  time,  p  is  pressure,  and  Re  =  UoL/v  is  the  Reynolds  number  based  on  a  characteristic 
length  L,  a  reference  velocity  U0,  and  the  kinematic  viscosity  v.  Equation  (1)  represents  the 
continuity  equation  and  Equation  (2)  represents  the  mean  momentum  equation.  The  equations  are 
written  in  tensor  notation  with  the  usual  summation  convention  assumed.  The  subscripts,  j  and  >7x, 
represent  the  covariant  derivatives. 

In  this  study,  the  two-layer  turbulence  model  of  Chen  and  Patel  (1988)  is  employed  to  provide 
closure  for  the  Reynolds  stress  tensor  u‘uJ  .  In  this  approach,  the  Reynolds  stresses  are  related  to  the 
corresponding  mean  rate  of  strain  through  an  isotropic  eddy  viscosity  vt: 


-u'uj  =  2v.Sij  —  —  giJk  (3) 

‘3 

with  V=|(g'*t/,i+gJWt);  k=LgW,  (4) 

where  SiJ  represent  the  contravariant  components  of  the  rate-of-strain  tensor  and  k  is  the  turbulent 
kinetic  energy.  Substituting  Equation  (3)  into  Equation  (2)  yields 


8U‘  ■■(  2  \ 

—  +  UJU‘+g‘J  p  +  -k 
dt  ’J  V  3  ) 


-2v,Sl'-jg',U'Jk=0. 


(5) 


The  quantity  1  /  Rt  =  1  /  Re  +  vt  represents  the  effective  viscosity.  In  the  fully  turbulent  flow  region 

away  from  the  solid  walls,  the  standard  k-s  model  is  employed  to  solve  the  transport  equations  for 
turbulent  kinetic  energy  k  and  its  dissipation  rate  s : 


dk  ■ 

—  +  UJk - 

dt  ’J 


f  l  ..  A 

\Kk  j 


-G+s=0 


—  +Uj£ 
dt  J 


f 


— Z  £ 


^  s  £2 

-Cel-G  +  C£2-  =  0, 

J,i  k  k 


(6) 

(V) 


where  the  eddy  viscosity  v,  and  the  production  term  G  are  given  by 

G=2v,gi„g„S»S™. 

£ 

The  effective  viscosities  in  Equations  (6)  and  (7)  are  taken  as  1 1  Rt  =  1 1  Re  +  v' /  a,s  and 
1  /  Re  =  1  /  Re  +  vt  /  ae ,  respectively. 


(8) 
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In  the  current  two-layer  approach  (Chen  and  Patel,  1988),  the  dissipation  rate  in  the  near  wall 
region  is  determined  from  the  turbulent  kinetic  energy  and  the  dissipation  length  scale  iE  to  account 
for  the  wall  effects: 

k3'2 

s  =  —  ;  l,=C,y[l-ew(-Ryl  A,)\.  (10) 

Using  this  relationship,  the  turbulent  kinetic  energy  can  be  determined  from  Equation  (6)  with  the 
following  eddy  viscosity  distribution: 

v,  =  c„^  (f  =  C,y[l-eM-KI 4)]-  <u> 

The  constants  C, ,  A  and  As  are  given  in  Chen  and  Patel  (1988)  are  chosen  to  yield  a  smooth 
distribution  of  eddy  viscosity  between  the  two  regions. 

The  above  equations  are  solved  numerically  using  the  finite-analytic  method  developed  by  Chen 
and  Chen  (1984),  Chen,  Patel  and  Ju  (1990),  Chen,  Bravo,  Chen  and  Xu  (1995),  and  Pontaza,  Chen 
and  Reddy  (2005).  A  detailed  description  of  the  finite-analytic  method  is  provided  in  Space  and 
Naval  Warfare  Systems  Center  Pacific  (SSC  Pacific)  Technical  Report  3293,  FANS-3D  Users  Guide 
Chen  and  Wang,  2016). 
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2.  SIMULATION  SCENARIOS  FOR  DDG  51  SHIP 


FANS  simulations  were  performed  for  a  USS  Arleigh  Burke  (DDG  51)  ship  as  shown  in  Figure  1 
under  two  different  water  depths  (10.0588  and  1 1 .5824  m)  and  two  different  propeller  rotating 
speeds  (26  and  5 1  rpms).  The  length  of  the  DDG  5 1  ship  is  142.04  m  (466  ft)  and  the  designed  draft 
is  9.4488  m.  The  diameter  of  the  twin-screw  propellers  is  5.4864  m  (18  ft),  and  the  center  of 
propeller  axis  is  located  at  5.7912  m  below  the  mean  water  level.  For  the  shallow  water  case  with 
10.0584  m  (33  ft)  water  depth,  the  underkeel  clearance  is  only  0.6096  m  (2  ft)  beneath  the  sonar 
dome  and  the  minimum  gap  between  the  propeller  tip  and  the  sea  bottom  is  1.524  m  (5  ft).  The 
propeller  rotating  speed  is  26  rpm  when  the  ship  speed  is  5  kts.  The  ship  speed  increases  to  10  kts 
when  the  propeller  is  rotating  at  5 1  rpm.  Detailed  information  of  the  DDG  5 1  ship  and  P4876 
propellers  are  summarized  in  Table  1. 


(b)  H  =  38  ft 


Figure  1.  DDG  51  ship  and  P4876  propeller  geometry. 
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Table  1.  Propeller  information  for  DDG  51  ship  in  FANS  simulation. 


Case  # 

1 

2 

3 

4 

Ship  length  L  (m) 

142.04 

142.04 

142.04 

142.04 

(466  ft) 

(466  ft) 

(466  ft) 

(466  ft) 

Ship  draft  (m) 

9.4488 

9.4488 

9.4488 

9.4488 

(31  ft) 

(31  ft) 

(31  ft) 

(31  ft) 

Water  depth,  H  (m) 

10.0584 

1 1 .5824 

10.0584 

1 1 .5824 

(33  ft) 

(38  ft) 

(33  ft) 

(38  ft) 

Underkeel  clearance  (m) 

0.6096 

0.6096 

2.1336 

2.1336 

(2  ft) 

(2  ft) 

(7  ft) 

(7  ft) 

Propeller  diameter,  D  (m) 

5.4864 

5.4864 

5.4864 

5.4864 

(18  ft) 

(18  ft) 

(18  ft) 

(18  ft) 

Distance  between  propellers  (m) 

9.8755 

9.8755 

9.8755 

9.8755 

(32.4  ft) 

(32.4  ft) 

(32.4  ft) 

(32.4  ft) 

Distance  from  ship  stern  to 

4.8768 

4.8768 

4.8768 

4.8768 

propeller  (m) 

(16  ft) 

(16  ft) 

(16  ft) 

(16  ft) 

Propeller  depth  (depth  of  the 

5.7912 

5.7912 

5.7912 

5.7912 

propeller  axis) 

(19  ft) 

(19  ft) 

(19  ft) 

(19  ft) 
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Table  1.  Propeller  information  for  DDG-51  ship  in  FANS  simulation.  (Continued) 


Case  # 

1 

2 

3 

4 

Distance  from  center  of  propeller 

4.2672 

4.2672 

5.7912 

5.7912 

axis  to  bottom 

(14  ft) 

(14  ft) 

(19  ft) 

(19  ft) 

Forward  Thrust  (N) 

47314 

47314 

175307 

175307 

Ship  speed  (knots) 

5 

5 

10 

10 

Propeller  rpm,  n 

26 

26 

51 

51 

Propeller  advance  coefficient 

1.082 

1.082 

1.103 

1.103 

Characteristic  velocity,  nD  (m/s) 

2.3774 

2.3774 

4.6634 

4.6634 

Reynolds  number  based  on 
propeller  diameter  D 

1.1 15  x  10" 

1.1 15  x  10" 

2.187  x  10' 

2.187  x  10' 

Reynolds  number  based  on  ship 
length  L 

3.262  x  10“ 

3.262  x  10“ 

6.245  x  10“ 

3.262  x  10“ 

2.1  FANS  MODEL  SIMULATION  RESULTS  FOR  DDG  51  SHIP 

Figure  2  shows  the  computational  domain  and  multi-block  overset  grids  used  in  this  study.  The 
overset  grid  system  consists  of  15  computational  blocks  and  7  phantom  grid  blocks  with  2,369,549 
grid  points  covering  half  of  the  solution  domain.  A  near-wall  spacing  of  5.4864  x  10"6  m  was  used 
near  the  sea  bottom  to  provide  accurate  resolution  of  the  turbulent  boundary  layer  flow.  Since  the 
first  grid  point  is  located  within  the  laminar  sublayer,  it  allows  us  to  calculate  the  shear  stresses  on 
the  seabed  directly  without  relying  on  the  wall-function  approximations. 

All  calculations  were  performed  for  4,000  time  steps  (i.e.,  100  propeller  revolutions)  using 
12  CPUs  on  a  Linux®  cluster.  For  the  5-kt  cases,  the  ship  travels  230.8  s  and  594  m.  When  the  ship 
speed  was  increased  to  10  kts,  it  took  approximately  1 17.7  secs  for  the  ship  to  travel  605  m  over 
100  propeller  revolutions.  The  simulation  results  clearly  indicated  that  the  propeller-induced  shear 
stresses  reached  a  periodic  pattern  in  less  than  50  propeller  revolutions. 


Figure  2.  Computational  domain  and  numerical  grids. 
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2.2  DISTURBED  VELOCITY  PROFILES 


Figure  3  shows  the  predicted  velocity  contours  and  velocity  vectors  adjacent  to  the  sea  bottom.  For 
completeness,  the  velocity  vectors  at  the  keel  plane  is  also  shown  in  Figure  4  to  provide  a  more 
detailed  understanding  of  the  three-dimensional  flow  field  induced  by  the  ship  motions.  For 
simplicity,  the  velocities  are  normalized  by  a  characteristic  velocity  V0  =  nD  given  in  Table  1,  where 
n  is  the  propeller  rotating  speed  (rps)  and  D  is  the  propeller  diameter  (m).  Note  that  the  ship  is 
traveling  in  the  negative  x-direction  on  an  earth-fixed  frame.  This  is  equivalent  to  a  positive  current 
in  the  x-direction  on  a  ship-fixed  reference  frame.  Figures  3  and  4  clearly  show  that  a  strong  flow 
acceleration  exists  beneath  the  sonar  dome  when  the  water  is  forced  to  pass  through  the  narrow 
underkeel  clearance  below  the  sonar  dome.  This  acceleration  resulted  in  positive  velocities  (in  the 
opposite  direction  of  the  ship  motion)  and  high  shear  stresses  beneath  the  sonar  dome. 

Figures  4(b)  and  4(d)  show  that  for  the  deep  water  cases  with  H  =  1 1.5824  m,  a  significant  portion 
of  the  flow  is  pushed  underneath  the  bow  due  to  local  flow  acceleration  around  the  sonar  dome.  This 
produces  a  large  flow  recirculation  region  (in  earth- fixed  reference  frame)  at  the  keel  plane,  with  a 
fairly  weak  return  flow  near  the  sea  bottom. 

When  the  water  depth  was  reduced  to  H  =  10.0584  m,  a  much  larger  resistance  pushes  the  flow 
beneath  the  sonar  dome.  Consequently,  most  of  the  surrounding  water  tends  to  move  laterally  around 
the  sonar  dome  and  the  flow  recirculation  near  the  seabed  was  confined  to  a  fairly  small  region 
immediately  downstream  of  the  bow,  as  shown  in  Figures  4(a)  and  4(c).  Figures  3(a)  and  3(c)  clearly 
show  that  the  ship  induced  a  strong  trailing  water  flow  (in  the  same  direction  of  the  ship  motion) 
beneath  the  ship  keel,  which  extends  beyond  the  propeller  plane  and  well  into  the  far  wake.  The  high 
velocity  (and  high  shear  stress)  regions  around  the  bow  and  mid-ship  are  induced  by  the  ship  hull 
movement,  but  not  directly  related  to  the  propeller  wash.  Note  that  the  effect  of  propeller  thrust  and 
torque  are  confined  to  the  ship  stem  and  wake  regions. 

Figure  5  shows  the  propeller-induced  velocity  distribution  along  the  center  plane  of  the  propeller 
axis.  The  axial  velocity  contours  at  five  selected  cross-sections  are  also  shown  in  Figure  6  to  provide 
a  better  understanding  of  the  swirling  flow  pattern  induced  by  the  propeller  rotation.  The  figure 
shows  that  the  propeller  rotation  induced  strong  swirling  flow  immediately  downstream  of  the 
propeller.  The  axial  flow  induced  by  the  propeller  thmst  force  remains  strong  for  more  than  15 
propeller  diameters  behind  the  ship  stem,  though  we  used  a  rather  coarse  grid  in  the  far  wake.  This 
strong  axial  flow  is  expected  to  carry  the  suspended  sediment  for  a  long  distance  downstream  of  the 
twin-screw  propellers. 


DDG-51,  V-5  kts ,  ft  PM =£6,  H-33  ft;  Time  Step-4000:  Time  -t  00 


DDS-51  V-6  kts  RPM-3&,  H-33  ji.  Time  Stef,- 4000  Time-  1 00 


(a)  Case  1:  V  =  5  kts,  H  =  10.0584  m,  co  =  26  rpm. 
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(b)  Case  2:  V  =  5  kts,  H  =  1 1 .5824  m,  co  =  26  rpm. 


DDC-5K  V- tQ  kts.  RPM- 5 I,  H-33  ft  Tim*  Step- 4000:  Time- 100 


(c)  Case  3:  V=10  kts,  H=1 0.0584  m,  co  =  51  rpm. 


(d)  Case  4:  V  =  10  kts,  H  =  11.5824  m,  co  =  51  rpm. 

Figure  3.  Longitudinal  velocity  contours  and  velocity  vectors  near  seabed. 
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(b)  V  =  5  kts,  RPM  =  26,  H  =  38  ft 


(a)  V  =  5  kts,  RPM  =  26,  H  =  33  ft 


(d)  V  =  10  kts,  RPM  =51,  H  =  38  ft 


Figure  4.  Velocity  profiles  at  the  keel  plane. 


(a)  Case  1:  V  =  5  kts,  H  =  10.0584  m. 


(b)  Case  2:  V  =  5  kts,  H  =  1 1 .5824  m. 


(c)  Case  3:  V  =  10  kts,  H  =  10.0584  m.  (d)  Case  4:  V  =  10  kts,  H  =  1 1 .5824  m. 

Figure  5.  Axial  velocity  contours  and  velocity  vector  plots  around  the  propeller. 
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In  addition  to  the  axial  flow  profiles  in  Figures  5  and  6,  the  swirling  flows  at  selected  cross- 
sections  are  also  shown  in  Figure  7  to  provide  a  complete  description  of  the  three-dimensional  flow 
field  induced  by  the  propeller  rotation.  The  figures  show  that  the  propeller-induced  swirling  flow 
patterns  are  quite  similar  immediately  downstream  of  the  propeller.  However,  the  propeller  swirl  is 
somewhat  stronger  in  the  far  wake  for  shallow-water  cases  with  H  =  10.0584  m.  Note  that  a  second 
pair  of  counter-rotating  vortices  are  near  the  center  plane  of  symmetry.  This  vortex  pair  was 
generated  in  the  narrow  gap  region  around  the  sonar  dome  and  remains  visible  in  the  far  wake. 

A  detailed  examination  of  the  velocity  profiles  near  sea  bottom  (see  Figures  3  and  4)  indicated  that 
the  propeller  wash  effect  is  negligible  for  the  deep  water  cases  since  the  swirling  flow  decreases 
quickly  in  the  radial  direction  away  from  the  propeller  tip  as  shown  in  Figure  7.  The  effect  of 
propeller  wash  grew  considerably  stronger  under  shallow- water  conditions  when  the  minimum  gap 
below  the  propeller  tip  was  reduced  from  3.048  m  (10  ft)  to  1 .524  m  (5  ft). 


DDG-51.  V-5  kts.  RPM-86.  H-3S  ft:  Time  Step-4000:  Time-tQO  DDG-St.  F=5  kts.  RPM=Z6,  H=38  ft:  Time  Step=4000.  Time  =100 


(a)  Case  1:  V  =  5  kts,  H  =  10.0584  m. 


(b)  Case  2:  V  =  5  kts,  H  =  1 1 .5824  m. 


Figure  6.  Propeller-induced  flow  field  at  selected  stations. 
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DDC  -  5 ; ,  V^Skts,  RPM=26,  H^33ft.  X/D^-65 


DDC-51.  V=5kts  RPM-26  H=38/t.  X/D^-65 


, ;  I  i  1 1  i !  1 1 U  iUuuaUiUUiUiUL 


DDC-51.  V=5kts.  RPM -26.  H=33ft,  X/D=-55 


U/(nD) 


DDC-51.  V=5kts.  RPM -26.  H=38ft,  X/D^-55 


(a)  Case  1 :  V  =  5  kts,  H  =  1 0.0584  m. 


(b)  Case  2:  V  =  5  kts,  H  =1 1 .5824  m. 
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(c)  Case  3:  V  =  10  kts,  H  =10.0584  m. 


(d)  Case  4:  V  =  1 0  kts,  H  =  1 1 .5824  m. 


Figure  7.  Propeller-induced  swirling  flows. 
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2.3  ESTIMATED  BOTTOM  SHEAR  STRESSES 


Figure  8  shows  the  shear  stress  distributions  on  the  sea  bottom  for  all  four  test  cases  considered  in 
the  current  study.  Note  that  different  color  bar  scales  were  used  since  the  bottom  shear  stresses  for 
Case  2  are  considerably  smaller  than  the  other  three  cases.  In  general,  the  bottom  shear  stress 
increases  with  the  propeller  revolutions  per  minute  and  ship  speed.  Under  deep-water  conditions,  the 
maximum  shear  stress  occurred  beneath  the  sonar  dome  due  to  strong  flow  acceleration  through  the 
narrow  passage  between  the  sonar  dome  and  sea  bottom.  High  shear  stress  regions  were  also 
observed  around  the  mid-ship  due  to  large  block  coefficient  of  the  DDG-5 1  hull  cross-section  area. 
The  propeller-induced  shear  stresses  are  not  as  high  as  those  induced  by  the  ship  motion. 
Furthermore,  the  propeller  wash  effects  are  confined  to  a  rather  small  region  directly  below  the  twin- 
screw  propellers.  For  shallow-water  cases  with  H  =  10.0584  m,  the  highest  shear  stress  also  occurred 
underneath  the  sonar  dome.  However,  the  shear  stresses  in  the  stem  region  are  also  very  high  due  to 
the  presence  of  strong  underkeel  current  (see  Figure  4)  induced  by  the  ship  motion. 


DDG-St.  V-S  lets.  PPM  =  26  H*33  ft:  Time  Step-3780:  Time-94. S 


DDC-St.  1 '=5  Us.  PPM  =  Sf  H-38  ft:  Time  Step- 3780:  Time =94. 5 


(a)  Case  1 :  V  =  5  kts,  H  =  1 0.0584  m.  (b)  Case  2:  V  =  5  kts,  H  =  1 1 .5824  m. 


DDG-St.  7m  to  Us.  RPM-Sl.  H-33  //.  Time  Stepm3780:  Time-94  S 


DDC-St.  V- tO  kts.  8PM  mS  I,  Hm38  ft.  Time  Step-3780:  Time-94.6 


(c)  Case  3,:V  =  10  kts,  H  =  10.0584  m. 


(d)  Case  4:  V  =  1 0  kts,  H  =  1 1 .5824  m. 


Figure  8.  Shear  stress  distribution  on  the  seabed. 


Figure  9  shows  the  surface  plots  of  seabed  shear  stress  distributions  under  different  water  depths 
and  different  ship  speeds.  For  clarity,  the  shear  stress  scales  were  adjusted  for  each  case  to  provide  a 
detailed  comparison  of  the  shear  stress  patterns  in  the  bow,  mid-ship,  propeller,  and  ship  stem 
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regions.  As  noted  earlier,  the  highest  shear  stress  occurred  beneath  the  sonar  dome  for  all  four  test 
cases  considered.  For  deeper  water  cases,  the  propeller  induced  shear  stresses  are  considerably 
smaller  than  those  induced  by  the  ship  hull  motion.  However,  the  bottom  shear  stress  distributions 
changed  drastically  when  the  water  depth  was  reduced  to  H  =  10.0584  m  with  a  very  small  underkeel 
clearance  of  0.6096  m.  The  high  shear  stresses  in  the  ship  wake  regions  were  induced  primarily  by 
the  trailing  water  in  the  narrow  gap  between  the  keel  and  seabed.  The  simulation  results  clearly 
demonstrated  that  the  blockage  effect  (i.e.,  block  coefficient)  of  the  ship  hull  in  shallow  water  is  the 
dominant  parameter  in  determining  the  sea  bottom  shear  stress  distributions. 


4 

3.5 


3- 


Y/D 


X/D 


(c)  Case  3,  V  =  10  kts,  H  =  10.0584  m 


(d)  Case  4,  V  =  1 0  kts,  H  =  1 1 .5824  m. 


Figure  9.  Surface  plots  of  shear  stresses  on  the  seabed. 


For  completeness,  enlarged  views  of  the  seabed  shear  stress  distributions  around  the  twin-screw 
propellers  and  the  ship  stem  regions  are  also  shown  in  Figure  10  to  provide  a  detailed  assessment  of 
the  propeller  wash  effects.  In  addition,  the  maximum  shear  stresses  at  various  locations  of  the  seabed 
were  also  summarized  in  Table  2.  The  propeller-induced  shear  stresses  are  much  smaller  than  the 
hull-induced  shear  stresses  under  shallow-water  conditions.  The  maximum  shear  stress  in  the  ship 
wake  exceeded  8.9  Pa  for  Case  3  when  the  ship  speed  was  10  kts.  Even  for  the  lower  speed  case  with 
V  =  5  kts  (Case  1),  the  bottom  shear  stress  in  Case  1  still  reached  nearly  2.5  Pa  in  the  ship  wake.  The 
maximum  shear  stresses  in  ship  wake  region  are  about  3  times  of  those  induced  by  the  propeller 
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rotation.  It  is  quite  obvious  that  the  large  blockage  effect  of  sonar  dome  under  shallow-water 
condition  is  the  primary  cause  of  the  strong  trailing  water  and  high  shear  stresses  in  the  narrow  gap 
between  the  ship  keel  and  sea  bottom.  However,  the  shear  stress  induced  by  the  ship  hull  is  strongly 
dependent  on  the  ship  size,  hull  form,  and  water  depth.  Therefore,  it  will  be  necessary  to  perform 
numerical  simulations  for  each  individual  ship  to  quantify  the  water  depth  effect  for  different  type  of 
ships. 


(c)  Case  3:  V  =  10  kts,  H  =  10.0584  m.  (d)  Case  4,  V  =  10  kts,  H  =  11.5824  m. 

Figure  10.  Surface  plots  of  seabed  shear  stresses  around  the  twin-screw  propellers. 


Table  2.  Maximum  shear  stresses  in  different  regions  of  the  seabed. 


Case  # 

1 

(Pa) 

2 

(Pa) 

3 

(Pa) 

4 

(Pa) 

Bow 

3.02 

1.09 

10.53 

3.57 

Mid-ship 

0.51 

0.39 

1.85 

1.46 

Propeller 

0.85 

0.18 

3.00 

0.85 

Ship  wake 

2.49 

0.14 

8.94 

0.37 
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3.  SIMULATION  SCENARIOS  FOR  TUGBOAT 
WITH  DUCTED  PROPELLERS 


FANS  simulations  were  also  performed  for  a  tugboat  with  two  ducted  propellers  as  shown  in 
Figure  1 1  at  a  constant  water  depth  of  9.144  m  (30  ft)  under  three  different  flow  conditions: 

(1)  propeller  blowing  to  open  water,  (2)  propeller  blowing  toward  a  pier  wall,  and  (3)  propeller 
blowing  parallel  to  a  pier  wall.  Since  the  detailed  geometries  of  the  tugboat  and  propellers  used  in  the 
experiments  are  not  available,  it  was  necessary  to  use  a  simple  barge-shaped  tugboat  as  shown  in 
Figure  10.  The  length  of  the  barge  is  28.65  m  (94  ft),  the  beam  is  10.36  m  (34  ft),  and  the  designed 
draft  is  3.35  m  (1 1  ft).  Also,  the  database  of  a  similar  ducted  propeller  was  used  to  represent  the 
actual  propellers  used  in  the  experiments.  The  propeller  diameter  is  scaled  to  2.286  m  (7.5  ft)  and  the 
hub  diameter  is  0.38  m  (1.25  ft).  The  outside  shroud  diameter  is  approximately  2.54  m  (100  in.).  The 
two  ducted  propellers  are  located  at  a  distance  of  15.24  m  (50  ft)  from  the  tugboat  stem,  and  at  a 
depth  of  4.88  m  (16  ft)  beneath  the  free  surface.  The  center-to-center  spacing  between  the  left  and 
right  propellers  is  4.88  m  (16  ft).  The  minimum  clearance  between  the  shroud  and  the  seabed  is 
approximately  3.00  m  (9.83  ft).  In  all  three  simulation  scenarios,  the  propellers  are  operated  at  200 
rpm  under  bollard-pull  condition  with  zero  forward  speed.  Detailed  information  of  the  barge-shaped 
tugboat,  ducted  propellers,  and  harbor  configurations  are  summarized  in  Table  3. 


Figure  11.  Tugboat  and  ducted  propeller  geometry. 
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Table  3.  Propeller  information  for  tugboat  in  FANS  simulation. 


Case  # 

1 

2 

3 

Ship  length  L  (m) 

28.65  (94  ft) 

28.65  (94  ft) 

28.65  (94  ft) 

Ship  beam  B  (m) 

10.36  (34  ft) 

10.36  (34  ft) 

10.36  (34  ft) 

Ship  draft  (m) 

3.353  (11  ft) 

3.353  (11  ft) 

3.353  (11  ft) 

Water  depth,  H  (m) 

9.144  (30  ft) 

9.144  (30  ft) 

9.144  (30  ft) 

Distance  from  ship  bow  to  pier 
wall  at  waterline  (m) 

1 .8288  (6  ft) 

Open  water 

1 .8288  (6  ft) 

Distance  from  ship  stern  to  pier 
wall  at  waterline  (m) 

Open  water 

30.48  m  (100  ft) 

Open  water 

Clearance  between  ship  sidewall 
and  pier  wall  (m) 

Open  water 

Open  water 

1.524  (5  ft) 

Underkeel  clearance  (m) 

2.997  (9.833  ft) 

2.997  (9.833  ft) 

2.997  (9.833  ft) 

Propeller  diameter,  D  (m) 

2.286  (7.5  ft) 

2.286  (7.5  ft) 

2.286  (7.5  ft) 

Distance  between  propellers  (m) 

4.8768  (16  ft) 

4.8768  (16  ft) 

4.8768  (16  ft) 

Distance  from  ship  stern  to 
propeller  (m) 

15.24  (50  ft) 

15.24  (50  ft) 

15.24  (50  ft) 

Propeller  depth  (depth  of  the 
propeller  axis) 

4.8768  (16  ft) 

4.8768  (16  ft) 

4.8768  (16  ft) 

Distance  from  center  of  propeller 
axis  to  bottom 

4.2672  (14  ft) 

4.2672  (14  ft) 

4.2672  (14  ft) 

Ship  speed  (knots) 

0 

0 

0 

Propeller  rpm,  n 

200 

200 

200 

Characteristic  time,  T0  (s) 

0.3 

0.3 

0.3 

Characteristic  velocity  U0  (m/s) 

1.016 

1.016 

1.106 

Reynolds  number  based  on 
characteristic  length  L0  (=  1  ft) 

2.647  x  10s 

2.647  x  10s 

2.647  x  103 

Reynolds  number  based  on 
propeller  diameter  D 

1.488  x  10' 

1.488  x  10' 

1.488x  10' 

We  considered  three  different  harbor  configurations  in  the  current  tugboat  propeller  wash 
study,  as  shown  in  Figure  12.  In  the  first  test  case,  the  tugboat  pushes  against  a  pier  wall  (with  a 
1.83-m  clearance  at  the  waterline)  and  the  ducted  propellers  blow  to  open  water.  For  the  second  test 
case,  the  propeller  blows  toward  a  pier  wall  located  at  30.48  m  (100  ft)  downstream  of  the  ship  stem. 
Note  that  the  distance  between  the  propeller  center  and  the  pier  wall  is  45.72  m  (150  ft)  since  the 
propeller  is  located  at  15.24  m  (50  ft)  upstream  of  the  ship  stem.  The  third  test  case  is  similar  to  Case 
1,  but  the  tugboat  is  aligned  in  parallel  to  another  pier  wall,  as  shown  in  Figure  12(c).  The  minimum 
distance  between  the  tugboat  sidewall  and  the  pier  wall  is  1.52  m  (5  ft).  As  noted  earlier,  the  rotating 
speed  (200  rpm)  of  the  ducted  propellers  and  the  water  depth  (30  ft)  remain  the  same  for  all  three  test 
cases. 
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(c)  Propeller  blowing  parallel  to  pier  wall. 

Figure  12.  Harbor  configurations  for  tugboat  propeller  wash  simulations. 
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3.1  FANS  MODEL  SIMULATION  RESULTS  FOR  TUGBOAT  WITH  DUCTED  PROPELLERS 


Figure  13  shows  the  computational  domain  and  multi-block  overset  grids  for  Case  3  with  propeller 
blowing  parallel  to  the  pier  wall.  The  overset  grid  system  consists  of  47  computational  blocks  and 
nine  phantom  grid  blocks,  with  7,070,832  grid  points  for  the  entire  solution  domain.  Fourteen  grid 
blocks  were  used  for  each  ducted  propeller  to  provide  detailed  resolution  of  the  propeller-induced 
flow  field  around  the  propeller  blades,  hub,  shroud,  and  near-wake  regions.  The  tugboat  is 
surrounded  by  a  single  body- fitted  grid  block,  and  the  far  field  is  divided  into  1 8  overlapping 
rectangular  grid  blocks.  The  computational  load  was  distributed  among  35  CPUs  on  a  Linux®  cluster. 


Figure  13.  Computational  domain  and  numerical  grids. 

For  simplicity,  we  performed  calculations  for  Cases  1  and  2  for  only  one-half  the  solution  domain 
using  Neumann  boundary  conditions  on  the  plane-of-symmetry  (y  =  0).  The  overset  grid  system  for 
Case  1  consists  of  24  computational  blocks  and  5  phantom  grid  blocks.  The  grids  around  the  ducted 
propeller  and  the  tugboat  are  identical  to  those  used  in  Case  3,  but  only  nine  computational  blocks 
are  needed  to  cover  the  far  field.  There  were  3,351,587  grid  nodes,  and  the  workload  was  distributed 
to  18  CPUs. 

The  propeller  and  tugboat  grids  for  Case  2  with  propeller  blowing  toward  a  pier  wall  are  also 
identical  to  those  used  in  Cases  1  and  3.  However,  we  had  to  refine  the  far- field  grid  in  front  of  the 
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pier  wall  to  provide  accurate  resolution  of  the  turbulent  boundary  layer  around  the  pier  wall.  The 
overset  grid  for  this  case  consists  of  29  computational  blocks  and  5  phantom  grid  blocks,  with  14 
rectangular  grid  blocks  covering  the  far  field.  We  used  4,736,735  grid  points  for  one-half  of  the 
solution  domain,  and  the  workload  was  distributed  among  23  CPUs  on  a  Linux®  cluster. 

For  all  three  test  cases,  a  near- wall  spacing  of  3.048  x  10"6  m  (10  5  ft)  was  used  next  to  the  sea 
bottom  to  provide  accurate  resolution  of  turbulent  boundary  layer  flow.  Since  the  first  grid  point  is 
located  within  the  laminar  sublayer,  it  allows  us  to  calculate  the  shear  stresses  on  the  seabed  directly 
without  relying  on  the  wall-function  approximations.  All  calculations  were  performed  for  12,500 
time  steps  (i.e.,  500  propeller  revolutions)  with  a  time  increment  of  0.012  sec. 

3.2  DISTURBED  VELOCITY  PROFILES 

Figures  14  and  15  show  the  predicted  velocity  contours  on  the  vertical  and  horizontal  propeller 
center  planes,  respectively,  for  Case  1  with  the  propeller  blowing  to  open  water.  For  completeness, 
the  velocity  contours  on  the  free  surface  are  also  shown  in  Figure  15  to  provide  a  more  detailed 
understanding  of  the  three-dimensional  flow  field  induced  by  the  ducted  propellers  under  a  bollard- 
pull  condition.  For  simplicity,  the  velocities  (U,  V,  W)  are  normalized  by  a  characteristic  velocity 
U0  =  nL0  given  in  Table  3,  where  n  is  the  propeller  rotating  speed  (rps)  and  the  characteristic  length 
L0  is  chosen  as  0.3048  m  (1  ft).  Also,  the  characteristic  time  is  defined  as  T0  =  L0/U0  =  1/n  such  that 
the  propeller  turns  one  revolution  over  one  characteristic  time.  For  current  simulations  with  the 
propellers  rotating  at  200  rpm,  the  corresponding  characteristic  time  and  velocity  scales  are 
T0  =  0.3  sec  and  U0  =  1 .016  m,  respectively. 

Figures  14  and  15  clearly  show  that  the  propeller  wake  extends  all  the  way  to  the  downstream  edge 
of  the  computational  domain  (about  76  m  behind  the  ducted  propellers)  in  less  than  300  propeller 
revolutions.  The  propeller  wake  flow  spreads  both  horizontally  and  vertically  over  several  propeller 
diameters.  Figure  14  shows  that  the  propeller-induced  velocities  are  considerably  stronger  near  the 
free  surface  due  to  local  flow  acceleration  in  the  narrow  clearance  between  the  propeller  blades  and 
the  tugboat  bottom  surface.  Note  that  the  tugboat  blockage  effect  is  strongly  affected  by  its  stem 
shape  and  cross-sectional  geometry  in  the  propeller  wake  regions.  The  blockage  effect  of  the  present 
barge-shaped  tugboat  tends  to  be  higher  than  the  other  ship-shaped  tugboats  since  the  block 
coefficient  of  a  rectangular  flat-bottom  barge  is  considerably  larger  than  other  tugboat  cross-sectional 
geometries. 

As  shown  in  Table  3,  the  minimum  distance  between  the  propeller  shroud  and  the  sea  bottom  is 
3.0  m,  which  is  only  about  1.3  times  of  the  propeller  diameter.  Due  to  the  relatively  small  underkeel 
clearance,  strong  propeller-induced  current  were  observed  near  the  sea  bottom  at  approximately  6  to 
10  propeller  diameters  downstream  of  the  propellers  during  the  initial  stage  of  the  propeller  wash 
simulation.  The  high  shear  stress  region  extends  to  more  than  20  propeller  diameters  downstream  at  a 
later  stage  as  the  propeller  wake  grows  considerably  longer  and  wider  after  t/To  =  400.  Note  that  the 
propeller-induced  flow  is  highly  three-dimensional  and  strongly  affected  by  the  tugboat  and  sea 
bottom  in  confined  harbor  with  shallow  water  depth. 
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Propeller  Blowing  To  Open  Water,  Time  Step  =  2500;  Time  =  1 00 


Propeller  Blowing  To  Open  Water ,  Time  Step  =5000;  Time=200 


Propeller  Blowing  To  Open  Water ,  Time  Step  =  7500;  Time=300 


Propeller  Blowing  To  Open  Water ,  Time  Step=  1 0000;  Time  =  400 


Prop  el  Lei'  Blowing  To  Open  Water ,  Time  Step=  12500;  Time  =  500 


Figure  14.  Axial  velocity  contours  on  the  vertical  propeller  center  plane  at  t/T0  =  100,  200,  300, 
400,  and  500  (T0  =  0.3  secs);  Case  1  with  propeller  blowing  to  open  water. 

In  the  second  test  case,  FANS  simulation  was  performed  for  the  same  tugboat  and  ducted 
propellers  in  Case  1  but  the  propellers  were  blowing  to  a  pier  wall  located  at  30.48  m  (100  ft) 
downstream  of  the  tugboat  stem.  Figures  16  shows  the  velocity  contours  on  the  vertical  propeller 
center  plane,  while  the  velocity  contours  on  the  horizontal  propeller  plane  and  the  free  surface  are 
shown  in  Figure  17.  Figure  16  shows  that  the  propeller- induced  flow  reached  the  pier  wall  at  about 
t/To  =  132,  and  was  deflected  by  the  pier  wall  and  spread  both  horizontally  and  vertically  along  the 
pier  wall.  This  produced  a  stagnation  flow  region  with  high  impact  pressure  on  the  pier  wall,  as 
shown  in  Figure  18. 
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Propeller  Blowing  To  Open  Water,  Time  Step=2S00.  Txme=lOO 
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(a)  Propeller  plane,  t  =  30  sec. 

Propeller  Blowing  To  Open  Water,  Time  Step-7500,  Time -300 
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Propeller  Blowing  To  Open  Water,  Time  Step-2500,  Time- 100 


(b)  free  surface,  t  =  30  sec. 

Propeller  Blowing  To  Open  Water ,  Time  Step-7500,  Time -300 


(c)  Propeller  plane,  t  =  90  sec. 


Propeller  Blowing  To  Open  Water,  Time  Step=  1 2500,  Time-500 
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(e)  Propeller  plane,  t  =  150  sec. 


(f)  free  surface,  t  =  1 50  sec. 


Figure  15.  Axial  velocity  contours  on  the  horizontal  propeller  center  plane  (left)  and  free  surface 
(right)  at  t/T0  =  100,  300,  and  500  (T0  =  0.3  sec);  Case  1  with  propeller  blowing  to  open  water. 


In  the  second  test  case,  FANS  simulation  was  performed  for  the  same  tugboat  and  ducted 
propellers  in  Case  1  but  the  propellers  were  blowing  to  a  pier  wall  located  at  30.48  m  (100  ft) 
downstream  of  the  tugboat  stem.  Figures  16  shows  the  velocity  contours  on  the  vertical  propeller 
center  plane,  while  the  velocity  contours  on  the  horizontal  propeller  plane  and  the  free  surface  are 
shown  in  Figure  17.  Figure  16  shows  that  the  propeller- induced  flow  reached  the  pier  wall  at  about 
t/T0  =  132,  and  was  deflected  by  the  pier  wall  and  spread  both  horizontally  and  vertically  along  the 
pier  wall.  This  produced  a  stagnation  flow  region  with  high  impact  pressure  on  the  pier  wall,  as 
shown  in  Figure  18. 
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Propeller  Blowing  Toward  Pier  Wall,  Time  Step=2500;  t/To=1 00 


Propeller  Blowing  Toward  Pier  Wall,  Time  Step=3300;  t/To=1 32 


Propeller  Blowing  Toward  Pier  Wall,  Time  Step=5000;  t/To  =  200 


Propeller  Blowing  Toward  Pier  Wall,  Time  Step=7500;  t/To=300 


Propeller  Blowing  Toward  Pier  Wall,  Time  Step-10000;  t/To-400 


Propeller  Blowing  Toward  Pier  Wall,  Time  Step=  1 2500;  t/To  =  500 


Figure  16.  Axial  velocity  contours  on  the  vertical  propeller  center  plane  at  t/T0  =  100,  132,  200,  300, 
400,  and  500  (T0  =  0.3  sec);  Case  2  with  propeller  blowing  toward  pier  wall 
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Propeller  Blowing  Toward  Pier  Wall.  Time  Step  =  2500.  t/T=100 


Propeller  Blowing  Toward  Pier  Wall.  Time  Step=2500 .  t/T=100 


(a)  Propeller  plane,  t  =  30  sec. 


Propeller  Blowing  Toward  Pier  Wall.  Time  Step  =  ?S00 :  t/T=30Q  Propeller  Blowing  Toward  Pier  Wall.  Time  Step  =  7500.  t/T=300 


(c)  Propeller  plane,  t  =  90  sec.  (d)  free  surface,  t  =  90  sec. 

Propeller  Blowing  Toward  Pier  Wall.  Time  Step=1 2500;  t/T=500  propeller  Blowing  Toward  Pier  Wall.  Time  Step=1 2500;  t/T-500 


U/Uo 


(e)  Propeller  plane,  t  =  150  sec.  (f)  free  surface,  t  =  150  sec. 

Figure  17.  Axial  velocity  contours  on  the  horizontal  propeller  center  plane  (left)  and  free  surface 
(right)  at  t/T0  =  100,  300,  and  500  (T0  =  0.3  sec);  Case  2  with  propeller  blowing  toward  pier  wall. 
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Propeller  Blowing  Toward  Pier  Wall,  Time  Step  =  4000;  t/To=  1 60 


Propeller  Blowing  Toward  Pier  Wall,  Time  Step =8000;  t/To  =  320 


Figure  18.  Pressure  contours  on  pier  wall  and  sea  bottom,  t/T0  =  160  and  320;  Case  2  with  propeller 
blowing  toward  a  pier  wall. 

FANS  simulations  were  also  performed  for  the  third  test  case  with  the  ducted  propellers  blowing 
parallel  to  a  pier  wall.  Due  to  the  asymmetric  harbor  configuration,  it  was  necessary  to  calculate  the 
flow  over  the  entire  solution  domain.  Furthermore,  unlike  the  counter-rotating  propellers  in  Cases  1 
and  2,  we  consider  two  identical  right-handed  propellers  here  with  both  propellers  rotating  in  the 
same  direction.  Figures  19  and  20  show  the  velocity  contours  on  the  vertical  propeller  planes  for  the 
right  (near  pier  wall)  and  left  (away  from  pier  wall)  propellers,  respectively.  Furthermore,  the 
velocity  contours  at  horizontal  propeller  center  plane  and  free  surface  are  also  shown  in  Figure  21  to 
provide  a  more  detailed  description  of  the  complex  three-dimensional  flow  induced  by  the  twin 
propellers.  It  is  clearly  seen  that  the  propeller  wake  flow  is  strongly  affected  by  the  presence  of  the 
parallel  pier  wall.  In  general,  the  velocity  induced  by  the  right  propeller  is  much  stronger  because  the 
propeller-induced  radial  flow  momentum  is  deflected  by  the  pier  wall  and  redirected  toward  the  axial 
direction  along  the  pier  wall.  On  the  other  hand,  the  wake  flow  induced  by  the  left  propeller  is 
somewhat  weaker  than  those  observed  in  Cases  1  and  2. 
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Propeller  Blowing  Parallel  to  Pier  Wall,  Time  Step=2500,  Time=lOO 


Propeller  Blowing  Parallel  to  Pier  Wall;  Time  Step=5000;  Time=200 


Propeller  Blowing  Parallel  to  Pier  Wall;  Time  Step  =  7500;  Time=300 


Propeller  Blowing  Parallel  to  Pier  Wall;  Time  Step=1 0000;  Time=400 


Propeller  Blowing  Parallel  to  Pier  Wall,  Time  Step=  12500;  Time=500 


Figure  19.  Axial  velocity  contours  behind  the  right  propeller  (near  pier  wall);  Case  3  with 
propeller  blowing  parallel  to  pier  wall. 


Propeller  Blowing  Parallel  to  Pier  Wall,  Time  Step=2500,  Time=1 00 


Figure  20.  Axial  velocity  contours  behind  the  left  propeller  (away  from  pier  wall);  Case  3  with 
propeller  blowing  parallel  to  pier  wall. 
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Propeller  Blowing  PoroUH  to  Pier  Wall;  Time  $iep=2$0O;  Ttme  =  1O0 
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(a)  Propeller  plane,  t  =  30  sec 


propeller  Blowing  Parallel  £$  Pier  Time  $tep=5O0B.  Time=200 
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(c)  Propeller  plane,  t  =  60  sec 


(b)  free  surface,  t  =  30  sec 

Propeller  Stowing  Pamllet  to  Pier  Boll;  Time  $tep=5O0Q.  Time=2QQ 


(d)  free  surface,  t  =  60  sec 


(e)  Propeller  plane,  t  =  90  sec 


propeller  B  louring  Parallel  to  Pier  if  oil;  Time  Step=f0000;  Time=40Q 


(f)  free  surface,  t  =  90  sec 


(g)  Propeller  plane,  t  =  120  sec 


(h)  free  surface,  t  =  120  sec 

Propeller  Blowing  Parallel  to  Pier  tftilL  Time  Step=t2SQ0.  Time =$00 


(j)  free  surface,  t  =  150  sec 


Figure  21.  Axial  velocity  contours  on  the  horizontal  propeller  center  plane  (left)  and  free  surface 
(right)  at  t/T0  =  100,  200,  300,  400,  and  500;  Case  3  with  propeller  blowing  parallel  to  pier  wall. 


For  completeness,  three-dimensional  views  of  the  propeller-induced  flow  at  the  horizontal 
propeller  center  plane  and  four  selected  axial  stations  are  also  shown  in  Figure  22  to  provide  a 
complete  description  of  the  propeller  wake  evolution.  The  figure  clearly  shows  that  the  swirling  flow 
behind  the  right  propeller  is  strongly  affected  by  the  pier  wall  as  well  as  the  left  propeller.  In  addition 
to  the  deflection  of  swirling  flow  momentum  by  the  pier  wall,  a  very  strong  interaction  also  exists 
between  the  left  and  right  propellers,  as  shown  in  Figure  23.  For  Cases  1  and  2,  the  propeller-induced 
swirling  flows  are  affected  primarily  by  the  flat  bottom  surface  of  the  tugboat.  For  the  co-rotating 
propellers  considered  in  Case  3,  however,  the  swirling  flow  momentums  in  the  overlap  region 
between  two  propeller  wakes  is  partially  suppressed.  This  resulted  in  a  deflection  of  the  weaker  left 
propeller  wake  (away  from  the  pier  wall)  toward  the  sea  bottom,  as  shown  in  Figure  22. 
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Propeller  Blowing  Parallel  to  Pier  Wall;  Time  Step=  1 0000;  t/To=400 


Propeller  Blowing  Parallel  to  Pier  Wall;  Tim.e  Step  =  1 0000;  t/To  =  400 


Figure  22.  Axial  velocity  contours  at  horizontal  propeller  center  plane  (top)  and  selected 
cross-sections  (bottom);  Case  3  with  propeller  blowing  parallel  to  pier  wall. 
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(a)  Case  1:  X  =  -12.19  m,  t  =  150  sec. 


(b)  Case  2:  X  =  -12.19  m,  t  =  150  sec. 

Propeller  Blowing  Parallel  to  Pier  WalL,  t/To  =  500 


Crossflow  Pr/po 


(c)  Case  3:  X  =  -12.19  m,  t  =  150  sec. 

Figure  23.  Propeller-induced  swirling  flows  (normalized  by  U0)  at  t/T0  =  500. 
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3.3  ESTIMATED  BOTTOM  SHEAR  STRESSES 


Figure  24  shows  the  shear  stress  distributions  on  the  sea  bottom  for  Case  1  with  propeller  blowing 
to  open  water.  In  general,  the  shear  stresses  are  high  during  the  initial  transient  period  (less  than  100 
propeller  revolutions)  after  impulsive  start  of  the  propeller  rotation.  As  the  propeller  wakes  grow 
longer  and  wider,  the  high  shear  stress  region  is  pushed  downstream  and  eventually  reaches  a  nearly 
periodic  pattern  after  400  to  450  propeller  revolutions.  At  t/T0  =  500,  the  high  shear  stress  region  is 
observed  between  5  to  10  propeller  diameters  downstream  of  the  propeller  and  the  maximum  shear 
stress  is  about  4.0  Pa  at  X/D  =  8.9.  Note  that  several  high-pressure  regions  exist  because  the  propeller 
wake  flow  is  highly  unsteady  in  confined  water  depth  under  bollard-pull  condition. 


Propeller  Blowing  To  Open  Water,  Time  Step=Z500,  Time=lOO  Propeller  Blowing  To  Open  Water,  Time  Step  =5000;  Txme=2QQ 


Propeller  Blowing  To  Open  Water,  Time  Step-7500;  Time -300 


Propeller  Blowing  To  Open  Water,  Time  Step=  10000  Time  =  400 


Propeller  Blowing  To  Open  Water,  Time  Slep-ttSOO;  Time  =  460  Propeller  Blowing  To  Open  if'ftSen  Time  Step- 12500;  Time  =  500 


Figure  24.  Shear  stress  distribution  on  the  seabed  at  t/T0  =  100,  200,  300,  400,  460,  and  500;  Case 
1  with  propeller  blowing  to  open  water. 
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Figure  25  shows  the  shear  stress  distributions  on  the  sea  bottom  for  Case  2  with  propeller  blowing 
to  a  pier  wall.  Before  the  propeller  wake  reaches  the  pier  wall,  the  shear  stress  pattern  for  Case  2  is 
similar  to  that  observed  in  Case  1  with  propeller  blowing  to  open  water.  After  the  propeller  wake 
impinges  on  the  pier  wall  at  about  t/T0  =  130,  another  high  shear  stress  region  was  developed  in  front 
of  the  pier  wall,  as  seen  in  Figure  25(c)-(f).  As  noted  earlier,  the  pier  wall  forced  the  propeller  wake 
to  spread  in  both  the  horizontal  and  vertical  directions.  Due  to  the  propeller  wake  impingement  and 
lateral  spread,  the  high  shear  stress  region  gradually  shifted  downstream.  At  t/T0  =  500,  the  high 
shear  stress  region  is  observed  near  the  pier  wall  with  a  maximum  shear  stress  of  about  4.2  Pa  at  X/D 
=  15.3  behind  the  duct  propellers.  Note  that  the  pier  wall  drastically  altered  the  shear  stress  pattern 
though  the  maximum  shear  stress  is  only  slightly  higher  than  that  observed  earlier  in  Case  1  with  the 
propellers  blowing  to  open  water. 


Propeller  Blowing  Toward  Pier  Wall,  Time  Slep*25Q0;  Time=-tOO  Propeller  Blowing  Toward  Pier  Wall.  Time  Step  =  5000 .  Time =200 


Propeller  Blowing  Toward  Pier  Wall,  Time  Step=7500 ;  Trme=300 


Propeller  Blowing  Toward  Pier  Wall „  Time  Step=  1 0000.  Time  =  400 


Propeller  Blowing  Toward  Pier  Wall.  Time  Step=1i500  Time  =  460 


Figure  25.  Shear  Stress  distribution  on  the  seabed  at  t/T0  =  100,  200,  300,  400,  460,  and  500;  Case 
2  with  propeller  blowing  to  pier  wall. 
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Figure  26  shows  the  shear  stress  distributions  on  the  sea  bottom  for  Case  3  with  propellers  blowing 
parallel  to  a  pier  wall.  Note  that  a  different  color  bar  scale  was  used  in  Figure  26  since  the  bottom 
shear  stresses  for  Case  3  are  several  times  higher  than  those  observed  in  Cases  1  and  2.  During  the 
initial  stage  of  simulation,  a  high  shear  stress  region  was  developed  around  the  tugboat  stem  region 
similar  to  those  observed  in  Cases  1  and  2.  As  the  propeller  wake  grew  longer  and  wider,  another 
high  shear  stress  region  was  developed  along  the  parallel  pier  wall,  as  shown  in  Figure  26(b)-(c). 


Propeller  Blowing  Parallel  to  Pier  Wall:  Time  Step*£500;  Time  =  1 00  Propeller  Blowing  Parallel  to  Pier  Wall;  Time  Step- 5000;  Time  =200 


Propeller  Blowing  Parallel  la  Pier  Wall;  Time  Siep=7500;  Time-m&OO  Propeller  Blowing  Parallel  to  Pier  Wall;  Time  Step  =  1 0000;  Time  =  400 


Propeller  Blowing  Parallel  lo  Pier  Wall;  Time  Step*  f  1 500:  Time**460  Propeller  Blowing  Parallel  to  Pier  Wall;  Time  Slepmt  2500;  Time=500 


Figure  26.  Shear  Stress  distribution  on  the  seabed  at  t/T0  =  100,  200,  300,  400,  460,  and  500; 

Case  3  with  propeller  blowing  parallel  to  pier  wall. 

As  shown  earlier  in  Figures  19  through  23,  a  very  strong  interaction  exists  between  the  left  and 
right  propellers  when  the  two  right-handed  propellers  are  rotating  in  the  same  direction.  Consequent¬ 
ly,  the  weaker  left  propeller  wake  (away  from  pier  wall)  was  pushed  closer  to  the  sea  bottom,  while 
the  stronger  right  propeller  wake  (near  pier  wall)  was  deflected  slightly  upward  to  the  free  surface. 
Note  that  the  weaker  left  propeller  wake  actually  produced  considerably  higher  shear  stresses  in  the 
near  field  because  it  is  much  closer  to  the  sea  bed  than  the  stronger  right  propeller  wake,  as  shown  in 
Figure  22.  In  the  far  field,  another  high  shear  stress  region  exists  adjacent  to  the  parallel  pier  wall.  A 
detailed  examination  of  the  propeller-induced  flow  field  indicates  that  the  second  high  shear  stress 
region  was  produced  by  the  stronger  right  propeller  wake.  At  t/T0  =  500,  the  high  shear  stress  regions 
is  observed  near  the  pier  wall  with  a  maximum  shear  stress  of  about  10.4  Pa  at  X/D  =  16.3  behind  the 
duct  propellers. 

A  detailed  comparison  of  the  velocity  and  shear  stress  distributions  for  Cases  1  through  3  clearly 
illustrates  that  the  sea  bed  shear  stresses  induced  by  the  ducted  propellers  are  strongly  affected  by  the 
harbor  configuration  as  well  as  the  propeller  rotating  directions.  Furthermore,  as  demonstrated  in  the 
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DDG  5 1  ship  and  P4876  propeller  wash  study,  the  shear  stresses  on  the  sea  bottom  also  depend 
strongly  on  the  water  depth,  ship  geometry,  underkeel  clearance,  ship  speed,  propeller  type,  and 
propeller  operating  conditions.  To  provide  detailed  shear  stress  distributions  for  different  ships  under 
site-specific  harbor  configurations,  it  is  desirable  to  perform  propeller  wash  study  by  solving  the 
Navier-Stokes  equations  directly  in  conjunction  with  advanced  turbulence  models. 
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